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We develop a general theory of electron-hole recombination in organic light emitting diodes that 
leads to formation of emissive singlet excitons and nonemissive triplet excitons. We briefly review 
other existing theories and show how our approach is substantively different from these theories. 
Using an exact time-dependent approach to the interchain/intermolecular charge-transfer within a 
long-range interacting model we find that, (i) the relative yield of the singlet exciton in polymers 
is considerably larger than the 25% predicted from statistical considerations, (ii) the singlet exciton 
yield increases with chain length in oligomers, and, (iii) in small molecules containing nitrogen 
heteroatoms, the relative yield of the singlet exciton is considerably smaller and may be even close 
to 25%. The above results are independent of whether or not the bond-charge repulsion, X±, is 
included in the interchain part of the Hamiltonian for the two-chain system. The larger (smaller) 
yield of the singlet (triplet) exciton in carbon-based long-chain polymers is a consequence of both 
its ionic (covalent) nature and smaller (larger) binding energy. In nitrogen containing monomers, 
wavefunctions are closer to the noninteracting limit, and this decreases (increases) the relative yield 
of the singlet (triplet) exciton. Our results are in qualitative agreement with electroluminescence 
experiments involving both molecular and polymeric light emitters. The time-dependent approach 
developed here for describing intermolecular charge-transfer processes is completely general and may 
be applied to many other such processes. 
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I. INTRODUCTION 



Charge recombination and photoinduccd charge- 
transfer lie at the heart of current attempts to con- 
struct viable optoelectronic devices using organic semi- 
conducting materials consisting of 7r-conjugated poly- 
mers or molecules. Charge recombination is the funda- 
mental process of interest in organic light emitting diodes 
(OLEDS). Electroluminescence (EL) in OLEDS results 
from, (a) the injection of electrons and holes into thin 
films containing the emissive material, (b) migration of 
these charges, which can involve both coherent motion 
on a single chain and interchain or intermolecular charge- 
transfer between neutral and charged species, (c) recom- 
bination of pJectrons and holes on the same polymer chain 
or moleculea'BEl. If the recombination leads to the singlet 
optical exciton, light emission can occur. If, on the other 
hand, the final product of the recombination is a triplet 
exciton, only nonradiative relaxation can occur in the 
absence of strong spin-orbit coupling. EL in OLEDS is 
of strong current rinterest, both because of applications 
in display devicesu'El and the potential for obtaining or- 
ganic solid state lasersa. The fundamental process that 
occurs in photoinduced charge-transfer is the exact re- 
verse of that in EL: optical excitation to the singlet ex- 
citon in a donor molecule is followed by charge separa- 
tion and migration of charge to a neighboring acceptor 
molecule. The latter process is of interest in photovoltaic 
applications! 



The fundamental electronic process of charge recombi- 
nation or separation is therefore of strong current inter- 
est. Especially in the context of EL in OLEDS, charge 
recombination has received both experimental and theo- 
retical attention (see below). The overall quantum effi- 
ciency of the EL depends on, (i) the fraction of the total 
number of injected carriers that end up as excitons on 
the same polymeric chain or molecule, (ii) the fraction 
of these excitons that are spin singlets, since only singlet 
excitons are emissive, and (iii) the fraction of singlet exci- 
tons that actually undergo radiative decay. In the present 
paper we focus on (ii), which determines the maximum 
possible EL efficiency. 

Formally, the charge recombination process can be 
written as, 

P + +P~ -^G + S/T (1) 

where P^ are charged polaronic states of the emissive 
molecule, G is the ground state of the neutral molecule, 
and S and T are singlet and triplet excited states of the 
neutral molecule. Eq. |l| indicates that both singlet and 
triplet excitons arc likely products of the charge recom- 
bination process. We shall denote the fraction of singlet 
excitons generated in OLEDS by the above recombina- 
tion process as r\. 

Early discussions of r\ were based on statistical argu- 
ments alone. Since electrons and holes are injected inde- 
pendently from the two electrodes, and since two spin-1/2 
particles can give three independent spin 1 states (with 
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Mg — -1, and +1) but only one spin state (Ms = 
0), it follows that rj is 0.25. Note, however, that this ar- 
gument is strictly valid only for noninteracting electrons, 
such that single-configuration molecular orbital descrip- 
tions of all eigenstates are valid. In such a case, the 
highest occupied and lowest unoccupied molecular or- 
bitals (HOMO and LUMO) are identical for the singlet 
and triplet excited states. Charge recombination (Eq. [I]) 
then involves merely the migration of an electron from 
the doubly (singly) occupied HOMO (LUMO) of P~ to 
the singly occupied (unoccupied) HOMO (LUMO) of the 
P + , for both singlet and triplet channels. The singlet 
channel and all three triplet channels of the charge re- 
combination process are equally likely within the MO 
scheme. If electrons are interacting, however, this simple 
single-configuration description breaks down, as all the 
states included in Eq. |l| are now superpositions of multi- 
ple configurations. There is no longer any fundamental 
reason for the singlet and triplet channels to be equally 
likely processes, and hence there is no reason for rj to be 
0.25. 

Experimentally, rj has been found to range from ~ 0.25 
- 0.6ffl§BEl in different materials. In OLEDS with the 
molecular species Aluminum tris (8-hydroxyquinolinc) 
(Alq3) as the emissive material Baldo et. al. have deter- 
mined -q ~ 0.22 ± 0.03, in .agreement with that expected 
from statistical arguments!! On the other hand, consid- 
erably larger rj ~ 0.45 has been found in derivatives of 
poly(parEuphenylenevinylene) (PPV) by Cao et al. and 
Ho et. al.BEJJ. Wohlgenannt et. al., using spin-dependent 
recombination spectroscopy, have determined the forma- 
tion cross-sections of singlet and triplet excitons, as and 
<jt , respectively, for a large number of polymeric materi- 
als (including nonemissive polymers in which the lowest 
two-photon state 2A g occurs below the optical 1B„ exci- 
ton), and found that (Js/ct is strongly material depen- 
dent and in all cases considerably larger than 1 (thereby 
implying that rj is material dependent apd much larger 
than 0.25)B. Jwlore recently, Wilson et alJl3 and Wohlge- 
nannt et. alB3 have shown that rj can depend strongly 
on the effective conjugation length, with values ranging 
from ~ 0.25 for small monomers to considerably larger 
than 0.25 for long chain oligomers. 

The oreti|Ca l hj ,. n has been investigated by a number of 
groupsEJiiJilj'tj'EZI including ourselves. There is general 
agreement that rj can be substantially greater than 0.25 
in 7r-conjugated polymers and that this is an electron 
correlation effect. There exist, however, substantial dif- 
ferences between the assumptions and formalisms that 
go into these theories. The goal of the present work is to 
develop a formalism that gives a clear physical picture of 
the electron-hole recombination and explains why rj sub- 
stantially larger than 0.25 is to be expected in organic 
polymeric systems. Ideally, since photoinduced charge- 
transfer is the exact reverse process of electron-hole re- 
combination, it should also be possible to extend our 
approach to photoinduced charge-transfer in the future. 
Brief presentation of our work has been made earlierllil, 



where, however, the emphasis was more on the experi- 
mental technique used by our experimental collaborators. 
Here we present the full theoretical details of our earlier 
work, provide a critique of the earlier theories and also 
report on the new and interesting results of our inves- 
tigation of external electric field effects on r), albeit for 
artificially large fields, and also on the role of nitrogen 
heteratoms in electron-hole recombination. Specifically, 
our theoretical approach involves a time-dependent for- 
malism, within which the initial state composed of two 
oppositely charged polarons is allowed to propagate in 
time under the influence of the complete Hamiltonian 
that includes both on-chain and interchain interactions. 
For the sake of completeness waafeo discuss other exist- 
ing theoretical approachesoEj£30, and their applicabil- 
ity to real systems. In particular, there exists a superficial 
similarity between the approach used in references [l6| [l?] 
and ours. For a physical understanding of the electron- 
hole recombination process it is essential (see below) that 
the difference between our approach and that used by the 
authors of reference [if], [l?] is precisely understood. 

The plan of the paper is as follows. In section ||, we 
present our theoretical models for intrachain and inter- 
chain interactions, and also discuss the model systems 
that are studied. In section III we present a brief cri- 
tique of the existing theories. A more extended discus- 
sion of the approach used by Shuai et. al.EjO is given 
in Appendix 1. In section IV we present the method of 
propagation of the initial state, while in section |v| we 
present our numerical results. In this section we also dis- 
cuss an alternate approach to the time propagation for 
the simplest case of two ethylenes that confirms the va- 
lidity of the more general approach used in section IV, 
and that also gives a physical picture of the recombina- 
tion process. While rj > 0.25 is found in our calculations 
with interacting electrons, the absolute yields of both sin- 
glet and triplet excitons are found to be extremely small 
with standard electron correlation parameters. We there- 
fore investigate the effects of the external electric field on 
these yields within a highly simplified model. It is found 
that for sufficiently large fields the yields with interact- 
ing electrons are as large as those with noninteracting 
electrons in the field-free case, and that in the relatively 
small field region rj continues to be greater than 0.25. In 
the very high field regime it is found that rj can even 
be smaller than 0.25. While the bare electric fields re- 
quired to see the reversal of the singlet-triplet ratio are 
rather large and therefore only of academic interest, if 
internal field effects are taken into account, it is pos- 
sible to envisage situations where the effective electric 
field is large enough to bring about such a reversal in 
the singlet-triplet ratio. Following the discussion of elec- 
tric field effects, we discuss how the cham-Length depen- 
dence of rj, as observed experimentallyOEj, can be un- 
derstood within our theory. We then consider the role of 
heteroatoms, especially in the context of molecular emit- 
ters. We show that in small systems with heteroatoms rj 
can approach the statistical limit, thus explaining quali- 
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tatively the monomer results oL Wilson et. al.li^, and the 
results of Baldo et. al. for Alq 3 El. The emphasis in all our 
calculations is on understanding the qualitative aspects 
of charge recombination and not on detailed quantitative 
aspects. Finally in section [v| we discuss the conclusions 
and scope of future work. 



II. THEORETICAL MODEL 

The goal of the present work is to provide bench- 
mark results for the charge recombination reaction which 
are valid for the strong Coulomb interactions that char- 
acterize 7r-conjugated systems. Accurate treatments of 
electron-electron interactions are not possible for long 
chain systems, and in this initial study we have therefore 
chosen pairs of short polyene chains, with 2-6 carbon 
atoms in each chain as our model systems. Since polyene 
eigenstates possess mirror-plane and inversion symme- 
tries, we shall henceforth refer to the ground state G 
(see Eq. |l|) as 1^4 g , and S and T as l 1 ^ and 1 3 B U , re- 
spectively. The model system containing two hexatrienes 
(12 carbon atoms overall) is the largest system that can 
be treated exactly at present within correlated electron 
models. 

Our approach suffers from two apparent disadvantages. 
First, polyenes and polyacetylenes are weakly emissive 
because the 2 1 A g state in these occur below the opti- 
cal l l B u state. This presents no problem as far as the 
analysis of the EL in emissive materials is concerned- 
as the spectroscopic technique of Wohlgenannt et. al.tiJ 
find a strong deviation of erg / <jt from 1 even in systems 
with energy ordering similar to that in polyenesEU (see 
results for PTV in this paper, for instance), and as we 
show in the following, this is a direct consequence of the 
large energy difference between the singlet exciton 
and the triplet 1 3 B U exciton, as well as the fundamental 
difference in their electronic structures. Both, in turn, 
are consequences of strong electron-electron interactions, 
which also characterize systems like PPV and poly- 
paraphenylene (PPP), as evidenced from the large differ- 
ence in energies between the singlet and triplet jeydtons 
in these systems, determined experimentallyEj'l!3oEj, as 
well as theoreticallyE3H. A second apparent disadvan- 
tage of our procedure is related to the limitation of our 
calculations to short systems. This prevents direct eval- 
uation of the chain length dependence of i]. We believe 
that this problem can be circumvented once the mech- 
anism of the physical process that leads to the differ- 
ence between singlet and triplet generation is precisely 
understood, and for this purpose it is essential that the 
electron correlation effects are investigated thoroughly 
using exactly solvable models. As we show later, our ap- 
proach gives a precise though qualitative explanation of 
the chain length dependence. 

Our model system consists of two polyene chains of 
equal lengths that lie directly on top of each other, sep- 
arated by 4 A. We consider the charge recombination 



process of Eq. Oj, and there are two possible initial states: 
(i) a specific chain (say chain 1) is positively charged, 
with the other (chain 2) having negative charge, a con- 
figuration that hereafter we denote as P 1 + P 2 ~, where the 
subscripts 1 and 2 are chain indices, or (ii) the super- 
position P^~P 2 ~ ± P2 P\ ■> i n the same notation. In our 
calculations we have chosen the first as the proper initial 
state, since experimentally in the OLEDS the symme- 
try between the chains is broken by the external electric 
field (we emphasize that the consequence of choosing the 
symmetric or antisymmetric superposition can be easily 
predicted from our all our numerical calculations that 
follow). Even with initial state (i), the final state can 
consist of both {I 1 A g ) x {1 1 B U ) 2 and (l 1 A g ) 2 (l 1 B u ) 1 in 
the singlet channel. The same is true in the triplet chan- 
nel, i.e., either of the two chains can be in the ground 
(excited) state. Hereafter we will write the initial states 
as \is) and |«t), where the subscripts S and T correspond 
to spin states S = and 1. We consider only the Ms = 
triplet state. The initial states are simply the product 
states with appropriate spin combinations, 

|is) = 2-V2(| P + )|p 2 - ) _ | P + )|p 2 -)) (2 ) 

\i T ) = 2- 1 / 2 (|P+)|P 2 : i ) + \P+ d )\P 2 -,)) (3) 

There exist of course two other initial triplet states 
with Ms = ± 1. The overall Hamiltonian for our com- 
posite two-chain system consists of an intrachain terms 
Hintm and interchain interactions Hi nte r ■ Additional in- 
teractions must be explicitly included to discuss exter- 
nal influences like the electric field etc. Hi ntra describ- 
ing individual .-chains is the Pariser-Parr-Pople (PPP) 
Hamiltonianc-aEa for 7r-electron systems, written as, 

Hintra = - } J t tj {a\ a a j(T + H.C.) + ^ e i n * + 

<ij>,<J i 

Uin h1 n hl + ^ v ij( n i - z i)( n j - z j) ( 4 ) 

i i>j 

where a\ creates a 7r-electron of spin a on carbon atom 

h ni.a — a \ a a i a is the number of electrons on atom i 
with spin a and rii — nj i<T is the total number of elec- 
trons on atom i, et is the site energy and Zi are the local 
chemical potentials. The hopping matrix element ty in 
the above are restricted to nearest neighbors and in prin- 
ciple can contain electron-phonon interactions, although 
a rigid bond approximation is used here. Ui and are 
the on-site and intrachain intersite Coulomb interactions. 

We use standard parameterizations for Hi ntra - The 
hopping integrals for single and double bonds are taken 
to be 2.232 eV and 2.568 eV, respectively and all the site 
energies of carbon atoms in a polymer with all equivalent 
sites are set to zero. We choose the Hubbard interaction 
parameter Uc for carbon to be 11.26 eV, and for the Vy 
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we choose the Ohno parameterizationB, 



ViA = 14.397 



28.794 

Ui + Ui 



(5) 



where the distance ry is in A, Vij is in eV and the local 
chemical potential zc for sp 2 carbon is one. It should be 
noted then when hetero atoms like nitrogen are present, 
the on-site correlation energy, the site energy and the 
local chemical potential could be different from those for 
carbon. For Hi nter , we choose the following form, 

H mter = -t± ^(4X:.<x + H.C.) + 

+X ± + <)(4x CT + H.C.) + 

V hi ( n * ~ z i) ( n 'j ~ z i' ) ( 6 ) 



In the above, primed and unprimed operators correspond 
to sites on different chains. Note that the interchain hop- 
ping t± is restricted to corresponding sites on the two 
chains, which are nearest interchain neighbors. The in- 
terchain Coulomb interaction Vi t j, however, includes in- 
teraction between any site on one chain with any other 
site on the other chain. In addition to the usual one- 
electron hopping that . o c cu rs within the zero differential 
overlap approximationoE3 we have also included a many- 
electron site charge-bond charge repulsion Xj_ (operating 
between nearest interchain neighbors only) that consists 
of multicenter Coulomb integrals. This term should also 
occur within Hi ntra , but is usually ignored there 
of its small magnitude, relative to all other terms 
In contrast, the t± in Hi nter is expected to be much 
smaller, and X± cannot be ignored in interchain pro- 
cesses, especially at large interchain separationsca. We 
have done calculations for both X± = and X± ^ 0. 



III. BRIEF CRITIQUE OF EXISTING 
THEORIES 



To put our work in the proper context we 
present a dis«iask»i|-e|f the existing theories of charge 
recombinationlljIialiaO in this section. The natures of 
Hintra within all these models are similar in the sense 
that they all incorporate intrachain Coulomb interac- 
tions, without which of course there cannot be any dif- 
ference between singlet and triplet generation. Follow- 
ing this, there is a fundamental difference between the 
models of references [l4|, |l^ on the one hand, and those 
of references [l6| |l7| and ours on the other. Within the 
theory of references [l4|, |l5|, there is no difference in sin- 
glet or triplet generation in the first stage of the charge- 
recombination process, which involves interchain charge- 
transfer. Within these models, interchain charge-transfer 
yields high energy singlet and triplet excited states of 



long chains that occur in the continuum, and the low- 
est singlet and triplet excitons result from relaxations 
of these high energy states. Differences in the relative 
yields of the lowest singlet and triplet excitons are con- 
sequences of differences in the intrachain relaxation pro- 
cesses in the singlet and triplet channels, that occur in the 
second stage of the overall process. In contrast, within 
our theoryO and the theory of references [l6] and [Tt], the 
lowest singlet and triplet excitons are generated directly 
from two oppositely charged polarons, and their differ- 
ent yields are consequences of the different cross-sections 
of the interchain charge-transfer reactions in the singlet 
and triplet channels. 

Within the model of Hong and MengQ, the contin- 
uum singlet state decays to the lowest singlet exciton, 
while the continuum triplet state decays to a high en- 
ergy triplet state T 2 consisting of a loosely bound triplet 
exciton, which then relaxes nonradiatively to the lowest 
tightly bound triplet exciton T%. The energy gap be- 
tween T 2 and T\ is large, and according to Hong and 
Meng, this nonradiative relaxation has to be a multi- 
phonon cascade process. The large energy gap and the 
multiphonon nature of the relaxation creates a "bottle- 
neck" in the T 2 — * T\ nonradiative transition, and spin- 
orbit coupling leads to intersystem crossing from T 2 to 
the singlet-exciton, thereby increasing the relative yield 
of singletstj. We believe that the key problem with this 
approach is that the model is in disagreement with what 
is known about the spectrum of triplet .states from triplet 
absorptions in 7r-conjugated-polymersE£l and theoretical 
solutions to the PPP model23. Experimentally, in PPA 
for instance, the lowest triplet occurs at about 1.55.eVE3, 
while in MEH-PPV this state occurs at - 1.3 eW3. The 
triplet absorption energy in these systems is about 1.4 
eV. Theoretically, the final state in triplet-absorption oc- 
curs slightly below the continuum bando, and this_is 
therefore the T 2 state (also referred to as the m 3 A g c3i). 
The energy region between T 2 and T\ (m 3 A g and 1 3 B U ) 
in the triplet subspace is not at all sparse, as assumed by 
Hong and Meng, but rather, within the correlated PPP 
Hamiltonian Hi ntra in Eq. ^, .fchis-|energy region contains 
numerous other triplet states£3'El Thus any nonradia- 
tive relaxation from T 2 to T\ in the realistic systems 
should involve a number of intermediate triplet states 
with small energy gaps between them, and therefore the 
phonon bottleneck simply will not occur. An additional 
problem with the model of Hong and Meng is that even in 
the singlet channel, generation of the lowest exciton from 
a continuum singlet state cannot be direct but can occur 
only through the m 1 A g loosely bound singlet excitonEl 
In principle, this can lead to a bottleneck even in the sin- 
glet channel. To summarize, we believe that the model 
of Hong and Meng is in disagreement with the known 
singlet and triplet energy spectra within the PPP model. 

Within the model of Kobrak and BittnerEl also po- 
laron pairs are formed on the single chain first. These 
authors take into account the electron-phonon interac- 
tions explicitly, and the two-particle states on a single 
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chain are allowed to evolve by interacting with a one- 
dimensional classical vibrational lattice. Different cross- 
sections for singlet and triplet excitons are found within 
the authors' model, and the difference originates from the 
difference in the mixing between the polaron and exciton 
states with different spin. The theory includes only the 
Coulomb interactions between the polaron charges and 
not the Coulomb interactions between all the electrons 
that appear in the PPP Hamiltonian. The theory also 
assumes large quantum efficiency for the generation of 
the high energy states with the two polaron charges on 
the same chain, starting from a state with the changes on 
different chains. A recent calculation by Ye et. al.EZl indi- 
cates very weak cross-sections for the generation of high 
energy 1 B U and 3 B U states starting from the initial state 
containingthe charges on different chains (see Fig. 8 in 
reference [17|) . This is supported also by our exact calcu- 
lations (see below). However, the calculations by Ye et 
al.l_J as well as ours are for relatively short chains, and 
further work is needed to test the validity of the model 
of Kobrak and Bittner. As we show in section [v], 77 > 
0.25 is predicted from considerations of the initial stage 
of interchain charge-transfer alone. Whether additional 
contributions can come from differences in the intrachain 
relaxation processes needs to be studied further. 

We now come to the work by Shuai et. al.E£lliZl, 
who, like us, have determined rj > 0.25 in oligomers of 
PPV from considerations of interchain charge-transfer. 
Precisely because of the apparent similarity of our ap- 
proaches, it is essential that we discuss the approach of 
Shuai et. al. in detail, since our ultimate goal is to ar- 
rive at a physical explanation of the greater yield of the 
singlet exciton than what is predicted from statistical 
considerations, and as we show later, the physical mech- 
anisms within references [16], |l7 and within our work are 
quite different. The quantity that is calculated in refer- 
ences |l^, ^ is erg/or, viz., the ratio of the formation 
cross-sections of the \ X B U singlet and 1 3 B U triplet exci- 
ton. For fast spin-lattice interaction, the expression for 
n in terms of 175 and ot can be written asO't^l, 



f] = °s/{os + 3ct t ) 



(7) 



and thus, for as /ax > 1, V > 0.25. 

Shuai et. al. consider the same Hi ntra as us, and 
Hinter that is similar (see below). The authors then use 
the Fermi "Golden Rule" approach to calculate erg and 
or- According to the authors, the cross-section ratio is 
given by, 

Osl°T = \{l S \H m ter\fs)\ 2 /WT\H mter \f T )\ 2 (8) 

where \is) and \it) are the singlet and triplet initial states 
(see Eqs. 2 and 3), and \fs) and |/t) are the corre- 
sponding final states, respectively. Since the interchain 
Coulomb interaction is diagonal in the space of the states 
considered in Eq. the authors ignore Vij in Eq. ^| but 
retain the other terms. Shuai et. al. find that for X± 
= in Eq. ^|, when the interchain charge-transfer is due 



to the hopping t± only, the right hand side of Eq. || is ~ 
1, a result we agree with (see Appendix 1). The authors 
then claim that for nonzero positive X± , and for positive 
ij_ (note negative sign in front of the one-electron term 
in Eq. ||), the right hand side of Eq. || can be substan- 
tially larger than 1. The authors calculated the matrix 
elements in Eq. || for pairs of PPV oligomers in parallel 
configuration using-approximate methods (singles-config- 
uration interaction!!^ and coupled-cluster methodO), and 
have found the right hand side of Eq. || to show diver- 
gent behavior over a broad range of X±/t± (see Fig. 1 
in reference [if] and Figs. 3, 4, 6 and 7 in reference [17]). 
Based on these calculations the authors conclude that a 
moderate to large Xj-is-pseniial for the experimentally 
observed large as / crTBt30t3'E3. 

This result is surprising, in view of the fact that the 
site charge-bond charge repulsion is spin-independent, ex- 
actly as the one-electron interchain hopping in Eq. |^. 
Since this question is intimately linked with the mecha- 
nism of charge recombination that we are after we have 
re-examined this issue by performing exact calculations 
for pairs of polyene chains with lengths N = 2, 4 and 
6. The conclusions from these exact calculations are de- 
scribed below. 

As discussed above, even with P^ P 2 as the ini- 
tial state (with, of course, appropriate spin func- 
tions) the final state contains two terms, with one 
of the two chains in the ground state and the other 
in the excited state. Instead of working with differ- 
ent superpositions of the final states we consider as 
to be proportional to \(is\H inte r\(l 1 A g ) 1 (l 1 B u ) 2 )\ 2 + 
\{is\H mter \{\ 1 A g ) 2 {l 1 B u ) 1 )\ 2 . Similarly, a T is taken 
to be proportional to \(iT\H in ter\(l 1 Ag) 1 (l 3 B u ) 2 )\ 2 + 
\{iT\H %nt er\{l 1 A g ) 2 {V i B u )i)\ 2 . As shown explicitly in 
the Appendix, the magnitudes of the matrix elements 
of the initial singlet [triplet] P^P 2 ~ with (l 1 A g )i(l 1 B u ) 2 
[(l 1 J 4 9 ) 1 (l 3 i? u ) 2 ] and (lM g ) 2 (l 1 B u ) 1 [(l 1 A 9 ) 2 (1 3 B U ) 1 ] 
are different for X± ^ 0, and hence the final states cannot 
be 1:1 superpositions of these configurations. Note that 
by taking the sums of the squares we exhaust all possibil- 
ities automatically. For the conclusions of references 
|i"7j to be valid the calculated as / or within Eq. |s| should 
now show strong dependence on Xj_/t± (as mentioned 
above divergent as /ax is implied in references |l6| , p~7| ) . 
Our exact results for the three different chain lengths are 
shown in Fig. 1 below, where we see that only for X±Jt± 
very close to 0.5 is as far, as calculated within Eq. is 
substantially different from 1. At all other X±/t± the 
RHS of Eq. H is very close to 1. Furthermore, except for 
X±/t± = 0.5 the chain length dependence of as far is 
weak. If we now recall that all chain length dependent 
quantities_(for example, optical and other energy gaps in 
polyeneso) exhibit strongest length dependence at the 
shortest lengths, the conclusion that emerges is that ex- 
cept for the unique point Xj_/t± = 0.5, as /ax remains 
~ 1 within the Golden Rule approach even in the long 
chain limit. 

In order to understand this difference from the results 
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of Shuai et. al.tS'tZl in further detail we present analytic 
results for the case of two ethylenes (N = 2) in Appendix 
1. These results are important in so far as they begin 
to give a physical picture for the charge recombination 
reaction, even as they indicate that the site charge-bond 
charge repulsion is not the origin of large rj. The analytic 
calculations also make the origin of the uniqueness of the 
point X±/t± = 0.5 absolutely clear. Indeed it is seen that 
precisely at this point both as and or, as defined in Eq. 
pL approach zero. More importantly, the chain length- 
independence, as suggested in Fig. 1 can be understood 
very clearly from the analytic calculations. Finally, it can 
also be seen from these calculations that had we taken 
the initial state to be the superposition P 1 + P 2 ~ i ^2 -Pi") 
instead of only one of these, the as/ar, as calculated 
from Eq. || would be exactly 1 for all X± /t± . 

Our basic conclusion then is that the Fermi Golden 
Rule approach is not valid for calculations of as /<jt or rj. 
This is to be expected also from a different consideration, 
viz., the Fermi Golden Rule approach is valid for calcula- 
tions of states that lie within a narrow band, whereas in 
the present case the energy difference between the initial 
and final states, and that between the singlet and triplet 
excitons are both much larger than t±_ and X±. The 
origin of the difference between our exact calculations 
of matrix elements and the approximate calculations of 
Shuai et. al. is harder to ascertain. One possibility is 
that the polaron wavefunctions are open shell, and ap- 
proximating these within mean field or limited CI could 
lead to wrong conclusions. 

In the following sections we therefore go beyond the 
Fermi Golden Rule approach to understand the origin of 
large rj. 



IV. TIME EVOLUTION OF THE POLARON 
PAIR STATE 

A straightforward numerical solution of 
Hintra + Hi nter will merely give the electronic structure 
of the composite two-chain system. Such a calculation 
does not contain any information about the relative 
yields of specific final states starting from the initial 
two-polaron states. Our approach therefore consists 
of propagating the initial state in time under the 
influence of the complete Hamiltonian, and monitoring 
the time-evolved state to obtain information about the 
final products. 

In principle, given a Hamiltonian, propagation of 
any initial state is easily achieved by solving the time- 
dependent Schrodinger equation. One could use the in- 
teraction picture to separate the nontrivial evolution of 
the initial state from the trivial component which occurs 
as a result of the evolution of the product oLihe eigen- 
states of the Hamiltonian of the subsystems 3 .]. In the 
context of the many-body PPP Hamiltonian such an ap- 
proach is difficult to implement numerically. This is be- 
cause the total number of eigenstates for the two-chain 



system is very large: the number of such states for two 
chains of six carbon atoms each is 853,776 in the M s 
— subspace. Obtaining all the eigenstates of the two- 
component system and expressing the matrix elements of 
Hi n ter in the basis of these eigenstates is therefore very 
intensive computationally. It is simpler to calculate the 
time evolution in the Schrodinger representation, deter- 
mine the time-evolved states, and project them on to the 
desired final eigenstates (for instance, \l 1 A g )i\l 1 B u ) 2 ). 
This is the approach we take. 

We first obtain the eigenstates \Py ), \P%) as well as 
the ptpduct states exactly in the valence bond (VB) 
basisE3 (in which the total spin 5 is a good quantum 
number) in order to avoid spin contamination. Following 
the time-evolution, however, we need to calculate over- 
laps of the time-evolved states with various final states 
(see below), which is cumbersome within the nonorthog- 
onal VB basis. After calculating the exact spin singlet 
and triplet initial states, we therefore expand these in an 
orthonormal basis that has only well defined total Ms 
value. 

Henceforth we refer to the initial states \ig) and 
collectively as ^(0) and the time-evolved states as ^(t). 
In principle, the time evolution can be done by operating 
on \& (0) with the time evolution operator, 

U(0,t)=exp(-iHt) (9) 

where H is the total Hamiltonian Hi ntra + Hi nte r- This 
approach would, however, require obtaining a matrix rep- 
resentation of the exponential time evolution operator, 
which in turn requires the determination of the pro- 
hibitively large number of eigenstates of the composite 
two-chain system. We can avoid this problem by using 
small discrete time intervals and expanding the exponen- 
tial operator in a Taylor series, and stopping at the linear 
term. Such an approach, however, has the undesirable 
effect of spoiling unitarity, and for long time evolutions 
would lead to loss of normalization of the evolved state. 
The way arou nd t his dilemma has been proposed and 
used by others 32 !! 33 ! in different contexts and involves us- 
ing the following truncated time-evolution scheme, 

(l + ifl-^*)tf(t + Af) = (l-iff^)*(t) (10) 

In the above equation, on the left hand side, we evolve 
the state at time (t + At) backwards by At/ 2 while on 
the right hand side, we evolve the state at time t forward 
by At/2. By forcing these two to be equal, we ensure 
unitarity in the time evolution of the state. It can be 
seen easily that this time evolution which is accurate to 
A|- is unitary. For a given many-body Hamiltonian and 
initial state, the right hand side of Eq. |l^ is a vector in 
the Hilbert space of the two-chain Hamiltonian. The left 
hand side corresponds to the action of a matrix on an 
as yet unknown vector, that is obtained by solving the 
above set of linear algebraic equations. Further details 
of the numerical procedure can be found in Appendix 2. 
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After each evolution step, the evolved state is projected 
onto the space of neutral product eigenstates of the two- 
chain system. The relative yield I mn (t) for a given prod- 
uct state \m, n) = \m)\\n)2 is then obtained from, 

J mn (t) = |(*(t)|m,n)| 2 (11) 

In our case the states \m, n) can be any of the final states 
of interest, viz., {(I 1 Ag)^! 1 B u ) 2 ) , |(l 1 A g ) 1 (l 3 B tl ) 2 ), etc. 
It is for efficient calculations of the overlaps (while at the 
same time maintaining spin purity) in Eq. [ll] that we 
expand our exact eigenstates of the neutral system in 
the VB basis to the total Ms basis. We emphasize that 
Imn(t) is a measure of the yield of the state \m, n) at 
time t and is not a cross-section. 



V. NUMERICAL RESULTS 

In this section we report the results of our calculations 
of recombination dynamics for for pairs of ethylenes, bu- 
tadienes and hexatrienes, both within the noninteracting 
Hiickel model (Z7j = Vy = X± = 0) and the interact- 
ing PPP model. Following this, we show the results of 
our investigation of electric field effects on the same sys- 
tems, discuss the chain length dependence of rj, and fi- 
nally present the numerical results for a model system 
containing nitrogen heteroatoms. The calculations for 
the noninteracting case provides a check of our numerical 
procedure, and the comparison between the noninteract- 
ing and the interacting model allows us to determine the 
effect of electron-electron interactions. 



A. Dynamics in the Hiickel Model 

While there is no difference in energy between singlets 
and triplets in the Hiickel Model, it is nevertheless pos- 
sible to have spin singlet and triplet initial states \is) 
and \it)i as well as singlet and triplet final states. In 
Fig. 2 we show the yield for the electron-hole recom- 
bination in the singlet channel, for pairs of ethylenes, 
butadienes and hexatrienes. The yields for the triplet 
channels are not shown separately in this case, - we have 
ascertained that these are identical to those in the sin- 
glet channel in this CHjSC, as expected. These calculations 
are for t±_ = 0.1 eV within Eq. ^. We note that the 
yields I mn (t) oscillate with time. This is to be expected 
within our purely electronic Hamiltonian, within which 
an electron or hole jumps back and forth between the 
two molecular species. These oscillations are the analogs 
of the Rabi oscillations 3 TP 5 ! that occur upon the stimu- 
lation of a system with light, where absorption of light 
can occur only with nonzero damping. Within our purely 
electronic Hamiltonian, complete transition to the final 
states can only occur in the presence of damping (for 
example, radiative and nonradiative relaxations of the fi- 
nal states), that has not been explicitly included in our 



Hamiltonian. The frequency of oscillation is higher for 
larger intermolecular transfer integral t±, as expected. 
The frequency of the oscillation also depends upon the 
size of the molecule and is lower for larger molecules (see 
below for an explanation of this). The equalities in the 
yields of the singlet and triplet excited states found nu- 
merically conforms to the simple free spin statistics which 
predicts that in the Ms = state formed from electron- 
hole recombination, the probability of singlet and triplet 
formation are equal. Since the Ms = ±1 cases always 
yield triplets, the spin statistics corresponding to 25% 
singlets and 75% triplets is vindicated in this case. 

Although the Hiickel calculations do not yield any new 
information, it is useful to pursue them further in order 
to arrive at a physical mechanism of the charge recom- 
bination process. To this end we have developed an al- 
ternate procedure for calculating the above dynamics for 
the smallest model system, viz., a pair of ethylenes. This 
alternate approach consists of expanding the initial state 
^(0) as a superposition of the eigenstates ipi of the com- 
posite two-chain system with eigenvalues £"i, 

|¥(0)>=X>hM0)) (12) 

i 

The evolution of the state ^(O) is now simply given by 
|tf (t)) = CillMO)} exp(-iEit/n) (13) 

i 

The yield I mn (t) in a given channel with final state \m 7 n) 
is then obtained from, 

(t) = \(m,n\*(t))\ 2 

= Y\d{m,n\^i (0)) eiq>(-iEit/h)\ 2 

i 

= £|c i (m,r# j (0))| a + 

i 

X^Z! 2 Re fe c i(ro,n|^(0))<^(0)|m,n>} 

i j>i 

x cos((K i - Ej)t/h) (14) 

The quantities (to, n\ipi(0)) are readily obtained from the 
eigenstates of the neutral one-chain subsystems and the 
composite eigenstates of the two-chain system. In Tables 
I and II we list the nonzero values of the coefficients Cj 
and the (to, n\ijji(0)} values for the case of two ethylenes. 
It is seen that sets of degenerate states of the compos- 
ite system together contribute equally to the singlet and 
triplet channels, although individual members of the set 
may contribute unequally. We have determined that the 
time evolution obtained from this approach is exactly the 
same as that obtained from the general method described 
in the previous section. 

The contribution arising from the right hand side of 
Eq. |l4] has been separated into time-independent and 
time-dependent parts. The latter comes about when- 
ever the two eigenstates in question are nondegenerate. 
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TABLE I: Significant c, =< *(0)|^i > and the < m,n\ipi > 
values and their product in the Hiickel model for a pair 
ethylenes in singlet channel. The index i corresponds to the 
index of 'significant' eigenstates of the total system and Ei 
the corresponding energy eigenvalue. 



i 


Ei (eV) 


Ci 


< m, n\ipi > 


< m,n\ipi > a 


2 


-4.3360 


0.3691 


0.1362 


0.0503 


3 


-4.3360 


-0.3373 


0.1138 


-0.0384 


4 


-4.1360 


-0.0171 


0.0120 


-0.0002 


5 


-4.1360 


-0.5000 


0.0058 


-0.0029 


6 


-4.1360 


0.4989 


0.0120 


0.0057 


7 


-4.1360 


-0.0285 


0.0059 


-0.0002 


8 


-3.9360 


-0.3558 


0.1266 


-0.0450 


9 


-3.9360 


0.3513 


0.1234 


0.0433 



TABLE II: Significant c, =< *(0)|V>i > and < m,n\ipi >, 
for the triplet channel, for a pair of ethylenes. 



i 


Ei (eV) 


< ¥(0)|V»4 > 


< m, n\tpi > 


< m,n |*i > Ci 


2 


-4.3360 


0.3373 


0.1138 


0.0384 


3 


-4.3360 


0.3691 


0.1362 


0.0503 


4 


-4.1360 


-0.0179 


0.0037 


-0.0001 


5 


-4.1360 


0.4985 


0.0180 


0.0090 


6 


-4.1360 


0.5005 


0.0047 


0.0024 


7 


-4.1360 


0.0251 


0.0170 


0.0004 


8 


-3.9360 


0.3513 


0.1234 


0.0433 


9 


-3.9360 


0.3558 


0.1266 


0.0450 



Furthermore, at t = the contribution from the time in- 
dependent part exactly cancels the contribution from the 
time dependent part. When the sign of the time depen- 
dent part becomes positive the two contributions add up 
to give the maximum yield of 0.25 in both the singlet and 
the triplet channels observed in the discrete calculations. 
The periodicity of the oscillation corresponds to the en- 
ergy difference between the two pairs of the degenerate 
states. This analysis could in principle be extended to 
the case of the larger systems but would be quite tedious 
in view of the larger Hilbert space dimensions. Note that 
the decrease of the oscillation frequency of I mn {t) with 
increasing chain length (Fig. 2) is explained within the 
above alternate procedure. The length dependence of the 
oscillation frequency originates from the smaller [Ei—Ej) 
in longer chains. 



B. Dynamics in the PPP model 

We now present our results for interacting electrons 
in Hi n t r a and Hi nter - In all cases for the interchain Vij 
we have chosen the Ohno parameters, and the interchain 
hopping t± =0.1 eV. For X±, we present the results of 
calculations with both X± = and 0.1 eV. In Figs. 3 (a) 



TABLE III: Significant a =< *(0)|^, > and the < m,n\ipi > 
values and their product for PPP model in the absence of 
electric field, for a pair of ethylenes in the singlet channel. The 
index i corresponds to the index of 'significant' eigenstates of 
the total system and Ei the corresponding energy eigenvalue. 





Ek 


Ci 


< m, n\tpi > 


Ci < m,n\ipi > 


4 


0.5295 


-0.0249 


-0.6992 


.0174 


5 


0.7328 


-0.0458 


-0.6953 


.0318 


11 


3.7748 


0.7066 


-0.0258 


-.0182 


13 


3.7844 


-0.7056 


0.0446 


-.0315 


29 


11.2503 


0.0082 


0.1020 


.0008 


30 


11.6379 


-0.0025 


0.1206 


-.0003 


32 


14.0483 


0.0050 


0.0028 


.00001 


34 


14.0611 


-0.0054 


0.0081 


-.00004 



and 3 (b) we show the plots of I mn (t) in the singlet and 
triplet channels for pairs of ethylenes, butadienes and 
hexatrienes, respectively, for the case of Xj_ = 0. The 
same results are shown in Figs 3 (c) and 3(d) for X± 
0. 1 eV. 

The most obvious difference from the Hiickel model is 
that the yields I mn (t) in both the singlet and triplet chan- 
nels are considerably reduced in the present cases. Two 
other points are to be noted. First, there is now substan- 
tial difference between the singlet and triplet channels, 
with the singlet yield higher in all cases. Second, the 
strong differences in singlet and triplet yields are true 
for both X± = and X± ^ 0.|— nOpis is in contradic- 
tion to the Golden Rule approachllaO, which ignores the 
energy difference between the l 1 ^ and the 1 3 B U . The 
only consequence of nonzero X±_ is the asymmetry be- 
tween the yields of {l 1 A g ) 1 (l 1 B u ) 2 and (l 1 A g ) 2 (l 1 B u ) 1 
in the singlet channels, and a similar asymmetry in the 
triplet channels. Further discussion of this asymmetry 
can be found in Appendix 1 . The overall conclusion that 
emerges from the results of Figs. 3 (a) - (d) is that 
nonzero electron-electron interactions substantially en- 
hances rj. 

In order to understand the above results in further de- 
tail we have also carried out the dynamics calculation for 
pairs of ethylenes according to Eq. [l4|. As in the Hiickel 
case these calculations yield the same results as the more 
general method. Our results for the wavefunctions of the 
composite two-chain system and the overlaps of the prod- 
uct eigenstates of the final neutral molecules with these 
are shown in Tables III and IV. The degeneracies in the 
eigenstates of the composite system that characterized 
the Hiickel model are now lifted, which is a known elec- 
tron correlation effect. What is more significant in the 
present case is that the composite state wavefunctions 
that have large overlaps with \I>(0) are now not the same 
ones that have large overlaps with the product wavefunc- 
tions of the final states. This is what reduces the yields of 
the charge-transfer processes in the PPP model, relative 
to the Hiickel model. 



9 



TABLE IV: Significant a =< *(0)|^» > and the < m, n\ipi > 
values and their product for PPP model in the absence of 
electric field, for a pair of ethylenes in the triplet channel. 





Ei 


c% 


<C TflTljUAlpi > 


d <T m Ti\ibi ~> 


2 


-2.7283 


-0.0215 


-0.6980 


.0150 


3 


-2.7238 


-0.0091 


0.6982 


-.0064 


10 


3.7697 


0.7068 


0.0060 


.0042 


12 


3.7775 


-0.7067 


0.0203 


-.0143 


19 


8.0804 


-0.0056 


0.1115 


-.0006 


20 


8.0875 


-0.0190 


-0.1114 


.0021 


31 


14.0475 


0.0051 


-0.0056 


-.00003 


33 


14.0515 


-0.0052 


-0.0023 


.00001 



Tables III and IV give a clear physical picture of the 
charge recombination process. For a large yield what 
appears to be essential is that the composite two- chain 
system must have at least some eigenstates which have 
simultaneously large overlaps with both the direct product 
of the initial polaronic states and the direct product of 
the pair of eigenstates of the neutral subsystems in the 
the chosen channel. This can be interpreted as a "tran- 
sition state theory" for the charge recombination reac- 
tion of Eq. [l]. Large overlaps with the initial polaronic 
pair states occur for the states 11 and 13 in the singlet 
channel (see Table III), and for the states 10 and 12 in 
the triplet channel (see Table IV). This is in contrast 
to the Hiickel case, where the large overlaps with the 
polaron pair wavefunctions were with the same compos- 
ite two-chain eigenstates. The overlaps of these specific 
two chain eigenstates are larger for products of singlet 
final states \l 1 A g )i\l 1 B u )2 than for triplet final states 
|l 1 A g )i|l 3 _B u )2, and this is what gives a larger yield for 
the singlet exciton. 

C. Effects of external electric field 

Our results in the previous subsection already indi- 
cate that ij can be substantially larger than 0.25 for the 
correlated electron Hamiltonian of Eq. ^. From compar- 
ison of Fig. 2 and Figs. 3 (a) - (d), we see however, 
that the relative yields I mn (t) are lower by orders of 
magnitude for interacting electrons. This is easily un- 
derstandable within time-independent second order per- 
turbation theory, within which the extent to which the 
initial polaron-pair state is modified is directly propor- 
tional to the matrix element of Hi nte r between the initial 
and final states, and inversely proportional to the zeroth 
order energy difference. Since the energy differences be- 
tween the polaron-pair states and the final neutral states 
are substantial within the PPP Hamiltonian, the yields 
are low. There are two possible interpretations of these 
results. First, the actual yields of excitons in OLEDS 
is indeed low, compared to the theoretical maximum for 
noninteracting electrons (recall that no direct compar- 



ison of the experimental light emission intensities with 
the theoretical maximum is possible). Second, the exper- 
imentally observed yields are influenced substantially by 
external factors ignored so far. We consider this second 
possibility here, and calculate within our time-dependent 
formalism the yields I mn (t) in the presence of an exter- 
nal electric field ( "external" in the following includes the 
effects of both the actual bias voltage as well as all in- 
ternal field effects). What follows may be thought of as 
overly simple, but nevertheless, we believe that it gives 
the correct physical picture. We first present our formal- 
ism and numerical calculations, and only then we discuss 
the interpretation of these results. 

As before, we consider pairs of molecules that are par- 
allel to each other, with the molecular chain-axes aligned 
parallel to the x-axis. The electric field is chosen along 
the y-axis, such that the total Hamiltonian now has an 
additional contribution, 

H field = Ej2((ni - l).Vi + {n'i - l).y$ (15) 

i 

In the above E is the strength of the electric field, and yt 
(y'i) gives the y-component of the location of the ith (i'th 
carbon atom in molecule 1 (2). We now perform our 
dynamical calculations with the complete Hamiltonian 
including H field- 

In Figs. 4 and 5 we show the effect of the external elec- 
tric field on the yield in the singlet and triplet channels 
for a pair of ethylenes. We see that in all the cases there 
is a strong nonlinear dependence of the yield on the ex- 
ternal field. In both the singlet and the triplet channels, 
we see sharp increases in the yields over a range of field 
strengths. The field strengths at which the increases in 
the yields occur are about two orders of magnitude larger 
than the experimental fields in the OLEDS, and we com- 
ment upon this below. Here we only observe that the 
field strength E over which the singlet yield is larger is 
smaller than field strength over which the triplet yield 
dominates. 

We have performed similar calculations for the longer 
chain systems, and in all cases the effects are the same, 
viz., there exists a range of field strength where a sud- 
den increase in the singlet yield occurs, while at still 
larger fields there occurs a similar jump in the triplet 
yield. In Figs 6 (a) and 6 (b) we have shown the singlet 
and triplet yields for field strengths of 0.3 V/A and 1.0 
V/A , respectively, for hexatriene. In general, for a given 
spin channel the threshold field strength decreases with 
the chain length (the threshold field for the singlet chan- 
nel decreases from 0.7 V/A to 0.3 V/A on going from 
ethylene to hexatriene, while the threshold field for the 
triplet channel decreases from 1.6 V/A to 1.0 VA ). The 
most important conclusions that emerge from these cal- 
culations are that, (a) macro scopically observable yields, 
comparable to the zero-field yields within the noninter- 
acting Hiickel model, are found for large fields, and (b) 
while the calculated n are greater than 0.25 for smaller 
fields, this is reversed with further increase in the field 
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TABLE V: In case of a pair of ethylenes the states with 
significant d =< 9(0)\t[>i > and < m,n\ipi >, for the PPP 
model with electric field of 0.7 V/A in the singlet channel. 



i 


Er (eV) 


a 


< m, n\ipi > 


a < m, n\tpi > 


4 


0.5153 


0.2086 


0.7629 


0.1591 


5 


0.6828 


0.3824 


0.4933 


0.1886 


11 


1.0479 


-0.9001 


0.3868 


-0.3482 


21 


6.5753 


0.0005 


-0.0055 


0.0000 


27 


11.2881 


-0.0016 


-0.0577 


0.0001 


30 


11.6946 


0.0040 


-0.1175 


-0.0005 



strength. 

In order to understand the origin of the increased yields 
over ranges of the electric field, we have analyzed the case 
of a pair of ethylenes extensively, within Eq. [l4|. Firstly 
it is worth noting that the geometry in which the field 
is introduced, the product states of the neutral Hamilto- 
nian are unaffected by the electric field. We also notice 
that the eigenvalues of the total Hamiltonian are not very 
sensitive to the external field. As in the field-free cases, 
we have obtained the projections of the eigenstates of the 
two-chain system on the initial state as well as the prod- 
uct of the final states, as a function of the applied electric 
field in both the singlet and the triplet channels. In Fig. 
7, we plot the coefficients (ipi\m,n) as function of the 
electric field for the singlet and the triplet channels. We 
see that there are several states that show strong vari- 
ation in both cases as a function of the field. However, 
when a product of these coefficients with < ^(O)]^ > is 
analyzed, the number of the states that simultaneously 
have a large value of these coefficients at the same elec- 
tric field is smaller. In Fig. 8, we plot the dominant 
coefficients of these projections, as a function of the ap- 
plied field. We note that only a few states have both 
projections simultaneously large. We also note that both 
the projections peak at the same field strength. It is this 
that leads to an abrupt increase in the yield at that field 
strength. 

The eigenstates of the full Hamiltonian that have large 
projections simultaneously to both the initial and the 
final states can in fact be expressed almost completely 
as a linear combination of the initial polaron product 
state and the final product state of the neutral system 
eigenstates. In Tables V and VI, we show the projec- 
tions of the eigenstates of the full Hamiltonian at the 
resonant electric field on to (i) the initial state and (ii) 
to the product of eigenstates of the neutral system for 
which resonance is observed. We note that there are a 
few eigenstates of the full Hamiltonian which have large 
coefficients for both projections. This seems to be in- 
dependent of the energy of the eigenstate of the total 
system. The energetics decide the period of oscillations 
and not the amplitude of the oscillations. 

We now come to our interpretations of the above nu- 
merical calculations. In all cases the applied fields in 



TABLE VI: In case of a pair of ethylenes the states with 
significant Ci —< »I/(0)|i/>i > and < m,n\ipi >, for the PPP 
model with electric field of 1.6 V/A in the triplet channel. 



i 


Ei (eV) 


a 


< m, n\tpi > 


Ci < m, Tl\lpi > 


2 


-2.8347 


0.5927 


0.3031 


0.1796 


3 


-2.7237 


0.0052 


0.9142 


0.0048 


5 


-2.5174 


0.8053 


-0.2179 


-0.1755 


20 


7.5928 


0.0089 


-0.0180 


-0.0002 


23 


8.0821 


0.0054 


0.1518 


0.0008 


24 


8.1387 


-0.0031 


0.0364 


-0.0001 


26 


10.1769 


0.0001 


0.0163 


0.0000 



our calculations are substantially larger than what is ex- 
pected from the externally applied voltage in OLEDS. 
Note, however, that our molecules are rather small, and 
the calculated threshold fields at which the effect be- 
comes observable decrease with the molecular size. In 
this context, it is worth recalling a previous exact cal- 
culation of electroabsorption for short finite polyenesEj. 
There the electric field was parallel to the chain axis (as 
opposed to perpendicular, as in the present case), and 
it was found that the calculated electroabsorption can 
simulate the experimentally observed behavior in long 
chain polymersEHH, provided the electric field used in 
the short chain calculation was larger by two orders of 
magnitude than the experimental field. This is because 
of the large energy gaps in short chains. We believe that 
a similar argument applies in the present calculations of 
interchain charge-transfer: the energy difference between 
the initial polaron-pair state and the final states is much 
larger in the small molecule-pair system than in the ex- 
perimental systems, even when oligomeric. The analogy 
to electroabsorption would then imply that the enhanced 
macroscopic yields would occur in the real systems at 
much smaller, perhaps even realistic fields. 

One final point concerns the geometry used in our cal- 
culations. In real OLEDS the relative orientations of 
the molecules of a given pair, as well as the orientation 
of the electric field with respect to individual members 
of the pair will both be different from that assumed in 
our simple calculations above. Electric fields that are 
nonorthogonal to the chain axis of a molecule will haw 
even stronger effects than found in our calculationsEa, 
while the random arrangements of the molecule pairs 
with respect to the field in the experimental systems im- 
plies that the range of field over which a given spin chan- 
nel dominates will be substantially larger than that found 
in our calculations. We therefore believe that a proper 
interpretation of our calculations is that in the experi- 
mental systems, there occur macroscopically large yields 
of both singlet and triplet excitons over a broad range 
of electric field. For small to moderate field strengths, 
the singlet channel dominates over the triplet channel. 
However, at still larger fields it is possible that this situ- 
ation reverses. Whether or not this higher regime of field 
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strength is experimentally accessible is a topic of future 
theoretical and experimental research. 



D. Chain length dependence 

We now discuss the chain length dependepcje-.of 77 as 
has recently been determined experimentallyli3ii3. From 
careful .measurements using different-techniques, Wilson 
et. al.E3 and Wohlgenannt et. alJij have established 
that 77 increases with conjugation length. Wilson et. al. 
have shown that while 77 is fifse to the statistically ex- 
pected 0.25 in the monomertU, it is substantially larger 
in the polymer. Wohlgenannt et. al. have shown that 
ct/cs increases linearly with the inverse of the conjuga- 
tion length. 

Within our numerical procedure, it is difficult to deter- 
mine the chain length dependence of 77 directly. This is 
because of multiple reasons, which include, (i) the limi- 
tation to rather small sizes, (ii) the necessity to integrate 
Imn(t) over a complete period in each case in Figs. 3 (a) 
- (d) to obtain the total yield over that period, and (hi) 
the difference in the periods for singlet and triplet chan- 
nels, as well as the differences among different singlet and 
triplet channels. We therefore present our discussion of 
chain length dependence within a simplified formalism 
that is consistent with our time-dependent procedure. 

Consider a transition (which could be of the charge- 
transfer type) between the states \k) and |s) of a two- 
state system, such that at time t — the system is in 
state \s) (c s (0) — 1, Cfc(0) = 0). We are interested in 
the yield |cfc(i)| 2 at a later time t due to a perturbation 
Vks- This is a standard textbook problemEa, and the 
time-dependent Schrodinger equation in this case is, 

ith— Ck = Vk s exp(iLu ks t) (16) 

For time-independent T4 S , as is true here (Vk s in the 
present case is simply Hi nter of Eq. ||) the above equation 
is easily integrated to give, 

In our case \s) — |P 1 + P 2 ~), with the appropriate spin 
combinations, and \k) — \(l 1 A g )i(\ 1 B u )2) for the sin- 
glet channel, and |fc) = \{\ 1 A g )\(l 3 B u )2) for the triplet 
channel (as usual, \k) can also have the chain indices 
reversed). We have already demonstrated (see section 
III and Appendix 1) that the matrix element (fc|F|s) is 
nearly the same for the singlet and triplet channels, ex- 
cept near the unique point X±/t± = 0.5. Ignoring the 
oscillation involving Uks we note that the relative yield of 
the singlet exciton is inversely proportional to the square 
of AE S = [E(P+)+E(P-)-E(l 1 A g )-E(l 1 B u )], while 
that of the triplet exciton is inversely proportional to 
the square of AEt, in which E^By) in the above is 
replaced with E(l 3 B u ) (note, however, that within the 



two-state approximation we have assumed that the sin- 
glet and triplet states that are of interest are the lowest 
singlet and triplet states; this can only be justified by the 
complete many-state calculations of the previous subsec- 
tions). We see immediately that this simple two-state 
formalism predicts higher singlet yield, since E^Bu) 
is considerably higher than E(l 3 B u ). Importantly, the 
chain- length dependence of 77 is also understood from the 
above. Both AEs and AEt decrease with increasing 
chain length. However, the ratio AEs / AEt also de- 
creases, because of the covalent character of the l 3 B u 
and the ionic character of the l 1 ^. This is seen most 
easily in the limit of the simple Hubbard model for the 
individual chains (zero intersite Coulomb interaction and 
zero bond alternation in Eq. ^), where AEs approaches 
and AEt approaches U in the long chain limit. 

We have calculated AEs and AEt exactly for all chain 
lengths N = 4 - 10 within the PPP-Ohno potential. 
While the ratio AEt/AEs shows the correct qualita- 
tive trend (viz., increasing AEt/AEs with increasing 
N) necessary for increasing 77 with increasing N, the ac- 
tual variation is small. This is to be expected, since our 
chain length variation is small, and the Ohno potential 
decays very slowly. With our limitation on N, it is nec- 
essary that the Coulomb potential is short range, such 
that we have the same Hamiltonian at all chain lengths, 
as is appro xima tely true for the experimental systems 
investigated!! 2 !!! 3 ]. We have therefore done exact calcula- 
tions of AEs and AEt for the extended Hubbard Hamil- 
tonian (Vij in Eq. |] limited to nearest neighbor interac- 
tion V) with parameters t%j — 1.08/jo and 0.92to for dou- 
ble and single bonds, V/to = 2 and U/to = 5 and 6. In 
Fig. 9 we show our calculated results for AEt/AEs for 
the two cases, for different N. In both cases, increasing 
AEt I AEs with increasing N indicates larger 77 for longer 
chain lengths. Energy convergences are faster with larger 
U, which explains the steeper behavior of AEt/AEs for 
larger U, and gives additional support to our argument. 



E. Role of heteroatoms 

The experiments by Baldo et. al.! and Wilson et. 
al.E 2 ] both indicate that in small molecular systems 77 
can be close to 0.25. This is in contrast to our results 
for ethylene (see Fig. 3). for which 77 is calculated to 
be substantially larger. One reason for this might be 
that the Coulomb correlation effects in thin film sam- 
ples are smaller than within the PPP Hamiltonian due 
to intermolecular interactions. The dominant effect, how- 
ever, is due to the heteroatoms in the molecules investi- 
gated by these authors, as we show below. Specifically, 
the site-energy (electronegativity) difference between the 
heteroatom and carbon atoms makes these systems closer 
to the Hiickel limit and this is what decreases 77. 

In order to compare with the model polyene systems we 
consider pairs of (CH=N)2 in the following calculations. 
The single chain Hamiltonian (Eq. 0) is then modified as 
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follows0. The Hubbard U for the nitrogen atoms, Un — 
12.34. eV. The local chemical potential zjy for nitrogen 
with 7r lone pair involved in conjugation is 2. Finally, 
nitrogen has site energy e = - 2.96 eV relative to that of 
the carbon atoms. There are two possible arrangements 
for the two chains in a parallel configuration, (i) carbon 
(nitrogen) on one chain lying directly above carbon (ni- 
trogen) atom on the other, and (ii) carbon atom on one 
chain lying above nitrogen atom on the second chain. We 
have chosen arrangement (i), - there is no fundamental 
reason for arrangement (ii) to have a very different rj. 

In Figs. 10 (a) and (b) we have plotted the I mn (t) 
for the singlet and triplet channels, respectively, for the 
case of X± = 0. Figs. 10 (c) and (d) show the same 
for X± = 0.1 eV. The most important conclusion that 
emerges from these calculations is that the relative yields 
of triplets are substantially larger in the present case, so 
much so that rj can be even close to the statistical limit 
of 0.25 (note that there are three triplet channels and 
the figures show the results for only one of these). We 
believe that these results give a qualitative explanation 
of the observation of Baldo et. al.H. Taken together with 
the chain length dependence of rj, as found in the pre- 
vious subsection, these results alpa explain qualitatively 
the observations by Wilson et. al.E3, since the same chain 
length dependence found in the case of simple polyenes 
should be true here also, although it is conceivable that 
rate of increase of rj with N here may be slower. 



VI. DISCUSSIONS AND CONCLUSIONS 

With a parallel arrangement of two polyene chains, 
we have shown that several experimentally observed 
qualitative features of the singlet-to-triplet yield ratios 
in 7r-conjugated systems can be understood within a 
well-defined total Hamiltonian for the two-chain system. 
While our model systems are rather simple, our theo- 
retical treatment of the charge-transfer process between 
the two chains is exact. We have given a full time- 
dependent approach to the interchain charge-transfer 
process, and have shown that in systems containing only 
carbon atoms, the overall yield of the singlet exciton is 
considerably larger than that of triplet excitons and 77 > 
0.25. This is a direct consequence of moderate electron- 
electron Coulomb interactions which has strong effects 
on both the energies and the wavefunctions of the singlet 
and triplet excitons. The mechanism of the exciton yields 
that emerges from our calculations is as follows. For 
large yields, it is essential that there exist excited states 
of the composite two-chain system whose wavefunctions 
have simultaneously large overlaps with the wavefunc- 
tion of the initial state consisting of polaron pairs, and 
the final state consisting of the two chains in the neutral 
states. Overlaps of the excited states of the two-chain 
system with final states in the singlet channel are con- 
siderably larger than for final states in the triplet chan- 
nel, and this is what gives a large yield for the singlet 



exciton. This is a consequence of the different natures 
of the singlet and triplet excitons, which are ionic and 
covalent, respectively, in the VB notation. Our result 
here is con^tanksath experiments on long oligomers and 
polymersB't2(lliHl3tj. Although our exact calculations are 
limited to short chains, within a two-state approxima- 
tion that is consistent with the full multi- level calculation 
we have shown that 77 increases with the chain length, 
in agreement with experimental observationaldllj. The 
two-state approximation gives an alternate explanation 
of the higher yield of the singlet exciton that is related 
to the singlet and triplet exciton binding energies, which 
are substantially different in 7r-conjugated polymers. Fi- 
nally, we have examined the role of heteroatoms, and 
have shown that in small molecular systems with nitro- 
gen as the heteroatom, rj is substantially smaller, and 
may be even close to the statistically expected value of 
0.25. The wavefunctions in this case, due to the strong 
electronegativity difference between the heteroatom and 
carbon atoms, are closer to the Hiickel limit, and this is 
what increases the relative yield of the triplet exciton. 
Our results, here successfully explain the difference-|be- 
tween Alq30 and heteroatom containing monomerstHI on 
the one hand, and polymeric systems on the other, and 
thereby provide additional strong support to our theo- 
retical approach. 

The time-dependent approach to the charge-transfer 
process developed here is completely general and can be 
applied to many other similar processes, for example, 
photoinduced charge-transfer, triplet-triplet collisions in 
OLEDS, etc. These and other applications are currently 
being investigated. Similarly, for a more complete un- 
derstanding of the chain length dependence of 77, we will 
investigate charge-transfer process within the density ma- 
trix renormalization group technique. 

While this manuscript was under preparation we re- 
ceived a manuscrip t (S. Karabunarliev and E.R. Bittner, 
cond-mat /020601 5 ) from E. Bittner that discusses the 
relative yields of singlet and triplet excitons within the 
context of intrachain processes (see section III) as op- 
posed to the interchain process discussed here. Although 
the approach of these authors is different from ours, they 
also find that the relative yields of singlet and triplet ex- 
citons are determined by their binding energies (smaller 
binding energies giving larger yields). It is not clear 
whether the approach used by these authors applies to 
molecule-based OLEDS. These authors have also investi- 
gated the effect of broken electron-hole symmetry, which 
is related to our calculations on chains of (CH=N)2- Our 
results here are different. While Karabunarliev and Bit- 
tner find even higher relative yield of singlet excitons 
(compared to electron-hole symmetric case) we find that 
•q here is smaller (see above). While a complete analy- 
sis of the electron-hole recombination must include both 
interchain and intrachain processes (and is a subject of 
future work in this area), we believe that this last result, 
when compared to experiments, justify our basic assump- 
tion that spin-dependence of the yields of excitons can be 
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understood largely within the context of intermolecular 
and interchain charge-transfer. 
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VIII. APPENDIX 1 

We present here detailed analytic calculations of the 
matrix elements of Hi nter for the case of two ethylenes. 
We believe that these calculations give clear understand- 
ings of the chain-length independence of the calculated 
as/err within the Fermi Golden Rule approach (Eq. ||) 
that was presented in section III (see Fig. 1). We also 
believe that even as these calculations show the inade- 
quacy of the Golden Rule approach they provide an indi- 
rect understanding of the actual mechanism behind large 
j] in long chain polymers. 

As in the rest of the paper we consider parallel ar- 
rangements of the ethylene molecules, with sites 1 and 2 
(3 and 4) corresponding to the lower (upper) molecule. 
Subscripts 1 and 2 that are assigned to wavefunctions 
describe the lower and upper molecule, respectively. The 
relevant single- molecule eigenstates, corresponding to the 
lower molecule then can be written as 



11%)! = (ci/V2)(ot t oI 4 + 4 T 4,i)|0) + 

(c 2 /V2)(<4 ;T 4 a - a| 4 a| 4 )|0) (16a) 

\1 X B U ) X = (l/V2)(aJ iT at lJr - 4 )T 4,|)I°) ( 16b ) 

ll 3 ^)! = (l/V2)(a t 1T 4 a + 4,;4, T )I > (16c) 

|P+) 1 = (l/V2)(a t 1)(7 +4 iCT )|0) (16d) 

|p-)i = (iA/2)(4 !T 4 4 4, ff -4,T4,i o L)i°> ( 16e ) 

In the above |0) is the vacuum for chain 1 and c\ and c 2 
are the coefficients of the ionic and covalent configura- 
tions in the lA g ground state that are to be determined 
by solving for the 2 x 2 A g subspace of ethylene within 
the PPP Hamiltonian (ci = c 2 = l/y/2 in the Hiickel 
Hamiltonian and the matrix elements in Eq. || in the sin- 
glet and triplet channels are exactly equal in this case). 
We have chosen the Ms = wavefunction for the 1 3 B U , 
but what follows is equally true for the Ms = ± 1 wave- 
functions. We have not assigned definite spin states to 
the charged polaronic wavefunctions, since the charged 
molecules can have either spin, and since different combi- 
nations of these spin states give the initial spin singlet or 



triplet product eigenstates for |P + P~). Note, however, 
the relative minus signs between the two configurations 
in |P~), as opposed to the relative plus signs between the 
two configurations in |P + ). This is what ensures that the 
product wavefunctions of the type \P^~ P 2 ~), with positive 
charge on molecule 1 and negative charge on molecule 2 
has odd parity with respect to the center of inversion on 
a single chain, and charge recombination can therefore 
only generate neutral states that have odd parity (for 
example, \l x A g ) 1 \l x B u ) 2 but not \l 1 A g ) 1 \2 1 A g ) 2 ). 

We consider the initial states \ig) and \ix) first, which 
are constructed from taking products of the polaronic 
wavefunctions given above. Since these product func- 
tions contain four terms each, and also since one of 
our goals is to arrive at a visual representation of the 
charge recombination process in configuration space, we 
have chosen not to write their explicit form but have 
given in Fig. 11 the wavefunctions in the VB notation, 
where a singlet bond between sites i and j is defined 

as 2~ 1 / 2 (4 f a tj.x — a l i a tj,T)|0), a triplet bond (with 
an arrow pointing from site i to site j) is defined as 
2~ 1 / 2 (4 | a tj.| + a \ x a tj.|)|0)i an d crosses correspond to 
doubly occupied sites 4f a tj JO). Given the initial and 
final states, it is now easily seen that Vij in -ffi«*™-4Eq. |^) 
plays no role within the Golden Rule approachll3tL3, since 
this term causes no transition between the initial and fi- 
nal states (note, however, Vij can play a significant role 
in the full dynamics calculation of section [V ) . The ma- 
trix elements of the remaining terms in H inter are now 
readily evaluated and these are given below, 

((lM 9 )i(l 1 B„) 2 |P mtei .|P 1 +P 2 -)s 
= -(ci/V2)(-i_L + 2X X ) - {c 2 /V2){-t ± + X ± ) (17a) 

((l 1 A 9 ) 2 (l 1 P„) 1 |P miel .|P+P 2 -) s 
= (ci/V2)(-t ± + 2X ± ) + (c 2 /V2)(-t± + 3X X ) (17b) 

((l 1 A g ) 1 (l 3 B u ) 2 \H mter \P+P^) T 
= -(ci/V2)(-tx + X ± ) - {c 2 /V2)(-t ± + 2X ± ) (17c) 

<(l 1 A s ) 1 (l 3 P u ) 2 |P mter |P+P 2 -) T 
= (ct/V2)(-t± + 3X X ) + (c 2 /V2){-t± + 2X±) (17d) 

Several points are to be noted now. First, for X± = 0, the 
squares of all the matrix elements are equal, and hence 
there is no difference between singlet and triplet genera- 
tion within the Golden Rule approach in this limit, and 
we agree on this with Shuai et. al£3 Second, however, 
defining overall as as the sum of the squares of the ma- 
trix elements in Eqs. 17(a) and and (b), and or as the 
sum of the squares of the matrix elements in Eqs. 17(c) 
and (d), respectively, we see that as/ur depends very 
weakly on X±/t± at all X±/t± except for X±/t± very 
close to 0.5, where — 1± + 2X± = and — 1± + X± and 
— 1± + 3X± have opposite signs. This is particularly so 
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for the calculated c\ and c 2 for PPP-Ohno parameters 
(ci = 0.5786, c 2 = 0.8156). We now examine the differ- 
ent terms in Eqs. 17(a) - (d) in detail. From Fig. 11 
we note that there are three classes of interchain elec- 
tron transfers: (i) charge-transfer between sites that are 
both singly occupied, leading to a doubly occupied site 
or an empty site (denoted by 1 + 1^2 + 0, where 
the numbers denote site occupancies) - or the exact re- 
verse process, (ii) charge-transfers of the type 1 + 0-^0 
+ 1, using the same notation, and (hi) charge-transfers 
of the type 2 + 1^1 + 2, again with the same no- 
tation. These three processes have different matrix ele- 
ments (-tj. + 2X±), (-i_L + X±) and (-t± + 3X±), re- 
spectively. The role of X± now becomes absolutely clear. 
Nonzero X± creates an asymmetry between the upper 
and lower molecule, leading to a difference between the 
yields of |l 1 A g )i|l 1 B tl ) 2 and \l 1 A g ) 2 \l 1 B u ) 1 , but it does 
not create a significant difference between as and ax- 

At exactly X±/tj_ = 0.5 terms containing (— t± + 2X±) 
in the matrix elements vanish, while the other terms are 
also small and of opposite signs. The singlet channel 
matrix elements now involve only c 2 . while the triplet 
channel matrix elements involve only c\. Since for re- 
pulsive Coulomb interactions c 2 > c\ , the sum of the the 
squares of the matrix elements here are larger for the sin- 
glet channel than for the triplet channel. This is what is 
reflected in our plot of Fig. ||. Note, however, that the 
calculated yields approach zero in both cases here. It is 
also clear from Eqs. 17 that this difference between the 
singlet and triplet channels persist over a very narrow 
region about X±/tj_ = 0.5. We therefore do not believe 
that this is of any relevance for realistic systems. 

Our final point concerns the chain-length independence 
of our results in Fig. || (except near Xj_/tj_ = 0.5). For 
arbitrary chain lengths there can occur only the three 
classes of interchain charge-transfers discussed above (1 
+ 1^2 + 0, 1 + 0^0 + 1 and 2 + 1^1 + 2). The 
detailed wavefunctions of longer chains are different, but 
the expectation values (nj^n^j) for the different wave- 
functions are nearly the same for fixed intrachain corre- 
lation parameters. Thus although in long chains there 
can in principle occur many more interchain hops that 
are of the type 1 + 1^2 + 0, such charge-transfers lead 
to additional double occupancies (relative to the over- 
all initial states) that are energetically costly because of 
electron correlation effects. Such charge-transfers there- 
fore make weak contributions to the overall interchain 
charge-transfer. The net consequence is the weak chain- 
length dependence found in Fig. 1 at all points other 
than X±/t± = 0.5. 

We believe that the above detailed calculation, aside 
from indicating the inapplicability of the Golden Rule, 
also indicates that the proper theoretical treatment must 
include the differences in the energies and wavefunc- 
tions of the final states, as indeed is done in our time- 
dependent calculations. 



IX. APPENDIX 2 

Details of the numerical procedure that were not dis- 
cussed in section |y| are given below. The charged as well 
as neutral eigenstates of Hi ntra for individual chains are 
obtained in the VB— representation by using a diagram- 
matic VB approaches with bit representation of the basis 
states. The eigenstates of a given spin S are obtained for 
Ms = S. The VB eigenstates are then transformed to the 
basis of Slater determinants with Ms — S by expanding 
the terms in each singlet pair and assigning an up-spin 
at each unpaired site with single occupancy. Thus, a 
triplet VB basis consisting of two singlet pairs and cor- 
responding to a function with Ms = 1 is expanded into 
four Slater determinants each with Ms = 1. To obtain 
eigenstate corresponding to other Ms values with Slater 
determinantal basis, we apply the S~ operator on the 
state, as many times as is necessary. 

We use the eigenstates of and \P~) to form the 
initial state of chosen spin in the form, 



*i,o(0) 
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2' + 2 
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(18) 
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2'~2 )l * l 2' + 2 ;2 



where the subscripts 1 and refer to the total spin of 
the initial state. The direct product of the states are 
expressed in the Slater determinantal basis of the com- 
posite system with the coefficient of the basis state 
in the composite system being given by 



<'k 



= did m < k\l*m > 



(19) 



l, in 



where di and d m are the coefficients of the basis states 
|Z) and \m) in the ground states of the subsystems 1 and 
2 respectively. The direct product itself is effected by 
shifting the 2n\ bits of the integer representing the basis 
state of system 1 with n\ sites to the immediate left of the 
2n 2 bits in the integer that represents the basis state of 
system 2 with n 2 sites. The resulting larger integer with 
2{n\ + 7i 2 ) bits correspond to an integer that represents 
one of the basis states of the composite system of («i+?i 2 ) 
sites. 

The evolution of the initial state involves solving the 
linear algebraic equations, Ax(t + At) = b, where, the 
matrix elements of the matrix A and the components of 
b are given by, 



Aij — 
b 



(dij + ifiHij-^-) 



At s 
~2 



At . 



(20) 
(21) 



The matrix A in the largest system we have studied is 
nearly of order one million and for reasonable conver- 
gence of the solution of the system of equations we need 
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a At of the order of Q.05eV/h which is typically 0.033 fem- 
toseconds, and this guarantees diagonal dominance of the 
matrix A. Thus, if one wishes to follow the dynamics for 
even as long as say 60 femtoseconds one needs to solve the 
linear system about 2000 times. This is rendered possible 
by the sparseness of the matrix A. For .efficient conver- 
gence, we use a small matrix algorithmcll which is very 
similar to the Davidson's algorithm for matrix eigenvalue 
problem. In the case of the largest system size, it takes 
about twelve hours on a DEC Alpha 333 MHz system to 
evolve the state by 60fs for a given channel. 

At the end of each iteration in the evolution of the 
state, we obtain the intensity or the yield in a pair of 



states in the neutral subsystem. The number of such 
pairs is enormous when we deal with say two systems of 
six sites each. The number of pairs in the singlet channel 
is 30,625 while that in the triplet channel is 33,075. We 
can reduce the pairs onto which we project the evolved 
state by restricting ourselves to a few low-energy states 
of the neutral subsystem. However, even in this case, 
the number of pairs could be rather large. To overcome 
this problem, we select only those pairs which have a 
minimum yield of say 10~ 4 at all times. This restriction 
when implemented judiciously leaves us with only a few 
pairs with significant yields. 
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Figure Captions 

Figure 1 : The ratio of the squares of the singlet and 
triplet matrix elements of Hi„t er (175 / ot according to Eq. 
||), as a function of X±/t± for pairs of ethylenes (circles), 
butadienes (squares) and hexatrienes (diamonds). 

Figure 2 : Yield in the singlet channel as a function 
of time, for pairs of ethylenes (top panel), butadienes 
(middle panel) and hexatrienes (bottom panel), within 
the simple Hiickel model (U = Vtj = X± = 0). Sig- 
nificant yield in all cases occur only for final states 
|(l 1 A+) 1 (l 1 B-) 2 > and |(l 1 B-) 1 (l 1 A+) 2 >, between 
which the yields are identical. Yields in the triplet 
channel \(l 1 A+) 1 (l 3 B+) 2 > and |(l 3 B+) 1 (l 1 yl+) 2 > are 
identical to those in the singlet channel. 

Figure 3 : Yields in the singlet and triplet chan- 
nels within the PPP Hamiltonian. In all cases the top 
panel corresponds to pair of ethylenes, the middle panel 
to pairs of butadienes, and the bottom panel to pairs of 
hexatrienes. (a) Singlet channel, t± = O.lcV, X± = 0; 
(b) Triplet channel, tj_ = O.lcV, X±_ = 0; (c) Singlet 
channel, t± = O.leV, X x = O.lcV; (d) Triplet channel 
t± = O.leV, X± = O.leV. Evolution in case of hexa- 
trienes is tracked for 20 fs while in other cases, the evo- 
lution is tracked for 60fs. Significant yields in singlet 
channel occurs only for final states \(l 1 A+)i(l 1 B~)2 > 
and |(l 1 B~)i(l 1 A+) 2 >, between which the yields are 
identical in (a) and (b) but different in (c) and (d). 
Similarly, yields in triplet channel are to the states 
|(liA+) 1 (l 3 S+) 2 >_andJ(l 3 B+) 1 (l 1 A+) 2 >, between 
which the yields are identical in (a) and (b) but different 
in (c) and (d). 

Figure 4 : Yields in the singlet channels (a) 
|(l%) 1 (lf) 2 ) ) (b) |(l 1 A 9 ) 2 (lf) 1 ), as a function of the 
electric field (V/A) and time (fs) . Here t± = 0.1 eV and 
X x = 0.1 eV. 

Figure 5 : Yields in the triplet channels (a) 
|(l 1 A 3 ) 1 (3f) 2 ), (b) |(l 1 A 9 ) 2 (3f) 1 ), as a function of the 
electric field (V/A) and time (fs). Parameters are same 



as in Fig. 4. 

Figure 6 : Yields in the singlet channel for pairs 
of hexatriene molecules, as a function of time (fs) with 
t± = O.leV and Xj_ — O.leV in an external electric field, 
(a) Singlet channel at 0.3 V/A, (b) singlet channel at 
0.42 V/Aand triplet channel at 1.0 V/A. 

Figure 7 : Evolution of significant < ipi\m, n > as 
a function of electric field (V/A), in case of the explicit 
time evolution of eigenvectors the PPP Hamiltonian for 
a pair of ethylenes in (a) singlet and (b) triplet channels. 

Figure 8 : < mn\ipi >< ^(0)|^i > plotted as a func- 
tion of electric field (V/A), for significant states T for 

(a) the singlet-singlet channel and (b) the singlet-triplet 
channel for a pair of ethylenes. The singlet-singlet chan- 
nel in (a) corresponds to \ms 1 > and \ns 2 > and the 
singlet-triplet channel in (b) corresponds to \ms t > and 
\n T >. 

Figure 9 : AEt / AE$ vs 1/N for the case of linear 
chains with "U-V" model Hamiltonian for the case of 
(i) U=5eV and V=2eV (squares) and (ii) U=6eV and 
V=2eV (circles). 

Figure 10 : Yields in the PPP model for the (CH=N) 2 
system, (a) singlet and (b) triplet channels with X± = 0; 

(c) and (d) singlet and triplet channels with X± = O.leV. 
The state to which the yield is significant in (a) is |SoSi > 
while in (b) it is to the state \SqT >. The yield to states 
I Si So > m singlet channel and |TSo > in triplet channel 
are identical to those for |SoSi > and \SqT > in (a) and 

(b) respectively. In (c) the yields to |SiSo > and |SoSi > 
are not the same and are shown separately. Similarly, in 

(d) yields to |TSo > and |SoT > are shown separately. 

Figure 11 : The initial (a) singlet and (b) triplet 
states \P1~P2) f° r the case of two ethylenes, and the 
result of operating with Hi nter - The upper (lower) two 
sites correspond to molecule 1 (molecule 2). The result 

(c) is a linear relationship between covalent triplet VB 
diagrams. 
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